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Abstract 

We consider the discrete time unitary dynamics given by a quantum walk on Z rf per- 
formed by a particle with internal degree of freedom, called coin state, according to the 
following iterated rule: a unitary update of the coin state takes place, followed by a shift on 
the lattice, conditioned on the coin state of the particle. We study the large time behavior 
of the quantum mechanical probability distribution of the position observable in Z d for 
random updates of the coin states of the following form. The random sequences of unitary 
updates are given by a site dependent function of a Markov chain in time, with the following 
properties: on each site, they share the same stationnary Markovian distribution and, for 
each fixed time, they form a deterministic periodic pattern on the lattice. 

We prove a Feynman-Kac formula to express the characteristic function of the averaged 

th 

distribution over the randomness at time n in terms of the n power of an operator M. By 
analyzing the spectrum of M, we show that this distribution posesses a drift proportional to 
the time and its centered counterpart displays a diffusive behavior with a diffusion matrix we 
compute. Moderate and large deviations principles are also proven to hold for the averaged 
distribution and the limit of the suitably rescaled corresponding characteristic function is 
shown to satisfy a diffusion equation. 

An example of random updates for which the analysis of the distribution can be per- 
formed without averaging is worked out. The random distribution displays a deterministic 
drift proportional to time and its centered counterpart gives rise to a random diffusion 
matrix whose law we compute. We complete the picture by presenting an uncorrelated 
example. 

1 Introduction 

Quantum walks are simple models of discrete time quantum evolution taking place on a 
(f-dimensional lattice whose implementation yields a unitary discrete dynamical system on 
a Hilbert space. The dynamics describes the motion of a quantum particle with internal 
degree of freedom on an infinite ci-dimensional lattice according to the following rules. The 
one-step motion consists in an update of the internal degree of freedom by means of a 
unitary transform in the relevant part of the Hilbert space, followed by a finite range shift 
on the lattice, conditioned on the internal degree of freedom of the particle. Due to their 
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similarity with classical random walks on a lattice, quantum walks constructed this way are 
often considered as their quantum analogs. In this context, the space of the internal degree 
of freedom is called coin space, the degree of freedom is the coin state and the unitary 
operators performing the update are coin matrices. 

Quantum walks have become quite popular in the quantum computing community in the 
recent years, due to the role they play in computer science, and in particular for quantum 
search algorithms. See for example [33], [4], [26], [28], [39], [5], [32] and in the review 
[36]. Also, quantum walks are used as effective dynamics of quantum systems in certain 
asymptotic regimes; see e. g. [14], [1], [33], [31], [11], [35], for a few models of this type, and 
[7], [10], [15], [17], [6] for their mathematical analysis. Moreover, quantum walk dynamics 
have been shown to describe experimental reality for systems of cold atoms trapped in 
suitably monitored optical lattices [24], and ions caught in monitored Paul traps [42]. 

The literature contains several variants of the quantum dynamics on a lattice as de- 
scribed above, which may include decoherence effects and/or more general graph, see e.g. 
the reviews and papers [4], [26], [3], [8]. In this work, we consider the case where the evo- 
lution of the walker is unitary, and where the underlying lattice is T, d with coin space of 
dimension 2d, which is, in a sense, the closest to the classical random walk. 

We are interested in the long time behavior of quantum mechanical expectation values 
of observables that are non-trivial on the lattice only, i.e. that do not depend on the 
internal degree of freedom of the quantum walker. Equivalently, this amounts to studying a 
family of random vectors X n on the lattice Z d , indexed by the discrete time variable, with 
probability laws P(X n = k) = Wk(n) defined by the prescriptions of quantum mechanics. 
The initial state of the quantum walker is described by a density matrix. 

As is well known, when the unitary update of the coin variable is performed at each time 
step by means of the same coin matrix, this leads to a ballistic behavior of the expectation 
of the position variable characterized by E w ^(X n ) ~ nV when n is large, for some vector 
V, and by fluctuations of the centered random variable X n — nV of order n, see e.g. [28]. 

The case where the coin matrices used to update the coin variable depend on the time 
step in a random fashion, a situation of temporal disorder, is dealt with in [21], see also 
[3]. All coin variables are updated simultaneously and in the same way, independently of 
the position on the lattice. This yields a random distribution W w {n), corresponding to 
the random variable X% which, once centered and averaged over the disorder, displays a 
diffusive behavior in the long time limit. 

If the coin matrices depend on the site of the lattice Z d but not on time, i.e. a case 
of spatial disorder, one expects dynamical localization, characterized by finite values of all 
moments, uniformly bounded in time n, and for (almost) all realizations. In dimension 
d = 1, this was proven in [20] for certain sets of random coin matrices, which were further 
generalized in [2]. See also [27], [38] for related aspects. The higher dimensional case is 
open. 

The situation addressed here is that of correlated spatio-temporal disorder. We consider 
random coin matrices which depend both on time and space in the following way: The 
random coin matrix at site x 6 Z, d and time n £ N is given by C%(x) = a x (io(n)), where 
{u)(j)}j<=fi is a temporal stationary Markov chain on a finite set of unitary matrices on C 2d , 
and 7L d 3 x — > a x is a given representation of Z d in terms of measure invariant bijections on 



2 



Q. In particular, ao =Id, the identity on Q, and r = {y € Z d s.t. a y = Id} forms a periodic 
sub-lattice of Z d . Therefore, at each site x G the sequence {C^(x)} je N is Markovian 
with a distribution independent of x, and at each time n G N, the set {C^(x), x € Z '} 
is T-periodic. This is a natural generalization of the case studied in [21] which displays a 
deterministic non trivial periodic structure in the spatial patterns of random coin matrices 
at each time step. 

This setup is an analog of the one addressed in [34], [22], [18], where the dynamics 
is generated by a quantum Hamiltonian with a time dependent potential generated by a 
random process. For quantum walks, the role of the random time dependent potential is 
played by the random coin operators whereas the role of the deterministic kinetic energy is 
played by the shift. 

We address the problem by an analysis of the large n behavior of the characteristic 
function of the distribution w.(n), $ n (y) = E, w ^(e iyXn ), where w.(n) = E(W w (n)) is the 
averaged quantum mechanical distribution on Z d , with initial condition po, a density matrix 
on Z 2 (Z d ) (g) C 2d . By adapting the strategy of [22], [18], inspired by [34], to our discrete time 
unitary setup, we first establish a Feynman-Kac type formula to express w. (n) in terms of 
(some matrix element) of the nP 1 power of a contraction operator M acting on an extended 
Hilbert space which involves a space of (density) matrices and the probability space of coin 
matrices. Then, we analyze the spectral properties of M, making use of the periodicity 
and invariance properties of o~ x which yield a fiber decomposition of a generalized Fourier 
transform of M. In turn, this allows us to provide a detailed description of the large n 
behavior of the characteristic function <3? n (y) in the diffusive regime y — > y/y/n, and at y 
fixed, in terms of the spectral data of M and their perturbative behavior. 

The foregoing is the main technical result of the paper, from which several consequences 
can be drawn, by arguments similar to those used in [21]. Under natural assumptions on the 
spectrum of M, the averaged distribution w.(n) displays a diffusive behavior characterized 
by the following data: a deterministic drift vector r 6 M. d and a diffusion matrix D, which 
we compute, such that, for n large and i,j = 1, 2, • • • , d, 

^w(n)i x n) - nr, E w{n) ((X n - nr)i(X n - nf)j) ~ n%. 

Moreover, we get convergence of the properly rescaled characteristic function of X n — nr, 
e -i[tn]ry/ v 7 "^^ (y j y/n), to the Fourier transform of superpositions of solutions to a diffusion 

equation of the form f Jd e~i {y \ n{p)y) dv / (2Ti) d with diffusion matrix D(p), p G T d , the d- 
dimensional torus. Also, we get moderate deviations results of the type 

F(X n -nr en (a+1 ^ 2 T) ~ e ~™ a mf.er A*M as n 0O> ( L1 ) 

for any set T £ M. d , any < a < 1, with some rate function A* : IR d — > [0, oo] we determine. 
Finally, we improve on [21] by establishing large deviations results for sets in a certain 
neighborhood of the origin, under stronger hypotheses. Informally, there exists an open 
ball B centered at the origin such that for all sets r G B n R d , 

¥(X n -nr enr)-e- ninfierI ' (l) as n -»• oo, (1.2) 

where A* : R d ->• [0, oo] is another rate function we determine. By Bryc's argument, [13], a 
central limit theorem for X n holds under the same conditions. 
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To complete the picture, we work out an example introduced in [21] where the distri- 
bution of coin matrices is supported on the set of unitary permutation matrices. This case 
allows us to analyze of the random distribution W w (n), without averaging over the disorder. 
We show that under our hypotheses, in this case W u (n) coincide with the distribution of 
a classical walk on the lattice, with increments being neither stationary, nor Markovian. 
Nevertheless, we can apply spectral methods as well to study the long times asymptotics of 
the corresponding random characteristic function, which allows us to get the the existence 
of a random diffusion matrix B w such that 

E w „ {n) ((X% - nr)i(X^ - nr),) ~ nD£, i, j = 1, 2, • • • , d, 

and whose matrix elements B^- are distributed according to the law of XfX^, where the 
vector X u is distributed according to JV(0, £), for a matrix E we determine. 

We also consider the completely decorrelated case where the coin matrices at each sites 
are i.i.d, i.e. a situation where no spatial structure is present in the pattern of coin matrices. 

Acknowledgements E.H. wishes to thank the CNRS and the Institut Fourier for 
support in the Fall of 2010, where this work was initiated, and A.J. wishes to thank the 
CRM for support in July 2011, where part of this work was done. 

2 General Setup 

Let H = C 2d (g> l 2 {Z d ) be the Hilbert space of the quantum walker in Z d with 2d in- 
ternal degrees of freedom. We denote the canonical basis of C 2d by {|r)} re/ d, where 
I± = {±1, ±2, . . . , ±d}, so that the orthogonal projectors on the basis vectors are noted 
Pt = \t)(t\, T £ I±- We shall denote the canonical basis of l 2 (Z d ) by {\x)} xeZ d, or by 
{S x } X £Z d - We shall write for a vector ip 6 tp = Yl,x&i d ^(^Ol^)' w here ip(x) = (x\tp) G C 2d 
and ^2 x€ %d ||V'( a ')llc 2d = ll^ll 2 ^ 00 • ^ e snan abuse notations by using the same symbols 
(•)•) for scalar products and corresponding "bra" and "ket" vectors on %, C 2d and l 2 (Z d ), 
the context allowing us to determine which spaces we are talking about. Also, we will often 
drop the subscript C 2d of the norm. 

A coin matrix acting on the internal degrees of freedom, or coin state, is a unitary matrix 
C G M2d{ < C) and a jump function is a function r : I± — > Z d . The shift S is defined on % by 

S = Pt ® \x + r(r)){x\. (2.1) 

xei d rei± 

By construction, a walker at site y with internal degree of freedom r represented by the 
vector |t) <8> \y) € % is just sent by S to one of the neighboring sites depending on r 
determined by the jump function r(r) 

S |r)® \y) = \T)®\y + r{T)). (2.2) 

The composition by C(y)<X>I, where the coin matrix C(y) is allowed to depend on the site y, 
reshuffles or updates the coin state so that the pieces of the wave function corresponding to 
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different internal states are shifted to different directions, depending on the internal state. 
The corresponding one step unitary evolution U of the walker on % = C 2d ®l 2 (li d ) is given 
by 

U = P T C(x)®\x + r{T))(x\. (2.3) 

x£l d r£J± 

Given a set of n > site-dependent unitary coin matrices Ck{x) G M2d(C), fc = 1, ■ ■ ■ , n 
and x G we construct an evolution operator U (n, 0) from time to time n, characterized 
at time k by Uk defined in (2.3) with {Ck(x)} xeZ d via 

U(n,0) = U n U n - 1 ---U 1 . (2.4) 

Let / : 7h d — > C and define the multiplication operator F : D(F) — > % on its do- 
main D(F) C % by {Ftjj)[x) = f(x)ij>(x), Vx G where V £ -C(-F) is equivalent to 

< 00 • Note that F acts trivially on the coin state. 
When / is real valued, F is self-adjoint and will be called a lattice observable. 

In particular, consider a walker characterized at time zero by the normalized vector ipQ = 
<fo <8> |0), i.e. which sits on site with coin state <po. The quantum mechanical expectation 
value of a lattice observable F at time n is given by (F),^ (n) = (ipo\U(n, Q)*FU(n, 0)?/>o)- 

A straightforward computation yields the following expression for the corresponding 
discrete evolution from time zero to time n 

Lemma 2.1 With the notations above, 

U(n,Q) = J2J2 J k( n )®\ x + k )( x l ( 2 - 5 ) 

x& d k& d 

where 

n—l n—2 

r k {n) = J2 Pr n C n {x+Y j r{T s ))p Tn _ 1 C n ^(x + Y j r(T s ))---P Tl C 1 {x) G M 2d (C) 

Ti,T2,---,T n £I± n s=l S=l 

(2-6) 

and J£(n) = 0, if Yl™=i r ( T s) / k. Moreover, for any lattice observable F, and any nor- 
malized vector tpo = y?o |0) ; 

{F)^(n) = <^|^*(n,0)F^(n,0)^)=X;/(fe)(^l^(nr^(n)^> 

= £/(fc)W fc (n), (2.7) 

where W k {n) = || J°(n)</? |lc2d safe/?/ 

£ W fc (n) = J] ||Jj?(n)^o||2 M = ||Y> |& = 1. (2.8) 
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Remark 2.2 We view the non-negative quantities {Wk(n)} n ^* as the probability distri- 
butions of a sequence of TL d -valued random variables {X n } ne ^* with 

Prob(X„ = k) = W k (n) = (ip \U{n,Oy(I®\k){k\)U(n,0)ip ) = || 4{n)<p \\ 2 c2d , (2.9) 

in keeping with (2.7). In particular, {F)^ (n) = E, Wk ^ n -j(f(X n )). We shall use freely both 
notations. 

Remark 2.3 All sums over k G Z fe are finite since J£(n) = i/max J= i i ^ \kj\ > pn, for 
some p > independent of x G since the jump functions have finite range. 

We are particularly interested in the long time behavior, n ^> 1, of (X 2 )^, (n), the expec- 
tation of the observable X 2 corresponding to the function f(x) = x 2 on Z d with initial 
condition ipo. Or, in other words, in the second moments of the distributions {Wk(n)} ne fq* . 

Let us proceed by expressing the probabilities Wk{n) in terms of the C^'s, k = 1, . . . , n. 
We need to introduce some more notations. Let I n (k) = {ti,--- ,r n }, where r\ G I±, 
I = 1, . . . , n and XT=i r ( r ') = ^- ^ n other words, /„(/c) denotes the set of paths that link 
the origin to k G Z d in n steps via the jump function r. Let us write (fo = ^2 re j ± a T \r). 

Lemma 2.4 

W fc (n) = ^ o^a^C^O) TlXnld^ro) x (2.10) 

T 0'{ T 1'"' . T n}€-ZnC0 

T-o.-f T i-'--' T n}e-rn(fc) 

s.t. t„=t4 



S-l S-l 

s=2 j=l j=l 



n s—L s—L 

II (<-i\<Z ( E Krj)) rj) (r s |c s ( £ rfo)) 



We approach the problem through the characteristic functions <3? ra of the probability 
distributions {W.(n)} n< =^* defined by the periodic function 

$ n (y) = E w(n) (e^) = W^nyy*, where y G [0,27r) d . (2.11) 



To emphasize the dependence in the initial state, we will sometimes write <l?n and/or 
W£°(n). All periodic functions will be viewed as functions defined on the torus, i.e. 
[0, 27r) d ~ T d . The asymptotic properties of the quantum walk emerge from the analy- 
sis of the limit in an appropriate sense as n — > oo of the characteristic function in the 
diffusive scaling 

lim $„(y/v^) (2.12) 

n^oo 

3 Correlated Markovian Random Framework 

We give here the hypotheses we make on the randomness of the model. 
Assumption C: 
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Let Q = {Ci,C2,--- ,Cf} be a finite set of unitary coin matrices on C 2d and let oj G 
fi N be a Markov chain with stationary initial distribution p and transition matrix P s.t. 
^(^j C) = Prob(a;(n + 1) = Cl w ( n ) = fi)i f° r an ^ G N. Let a be a representation of Z d , 
x i-> ctj;, in terms of measure preserving maps : — >■ 17 such that p(a x () = p(C) an d 

P( ff ,C ) ^) = P(C, J? ). 

Remarks 3.1 ^ This is equivalent to saying that the paths of a x (u(-)) have the same 
distribution as the paths ofcj(-), for all x £ Z d . 

ii) Because x ^ a x is a representation ofZ d , a x is a bisection over the finite set fi for any 
x G Z d and o~q = Id. Moreover, the finite set of bijections {o~ x } xe %d must commute with one 
another. 

Hi) Let r = {x G Z d s.t. a x = Id}. Then a x = o~ y is equivalent to x — y G T. If 
g G N* denotes the cardinal of the group {o^l^es*? then for any j G {1, • • • , d}, the vector 
(0, • • • , 0, g, 0, ■ ■ ■ , 0) T G Z d , where g sits at the j 'th slot, belongs to V. Hence the lattice Y 
is of dimension d. 

iv) We choose Br C Z d such that G .Br and o~\b v is a bisection on the set of bijections of 
Q. For any x G Z d , we have a unique decomposition x = xo + 1], with xq G -Br and r\ G V. 

We consider the random evolution obtained from sequences of coin matrices defined on 
site x G Z d at time n > by 

C£(s) = a x (u(n)). (3.1) 

This means that while the coin matrices at different sites have all the same distribution 
as C r ^(0) = oj(n), they can take different correlated values depending on a x . 

It is more natural in this setting to carry out the analysis in terms of density matrices. 
The set of density matrices, T>M, consists in trace one non negative operators on C 2d <8> 
l 2 (Z d ). Any bounded operators on U = C 2d (g> l 2 {Z d ) can be represented by its kernel as 

p= P( x iV)®\ x ){yl where /9(x,y) G M 2d (C). (3.2) 

A non-negative operator p on % is trace class iff 

\\p(x,x)\\<oo. (3.3) 

x ez d 

We say that p belongs to l 2 (Z d x Z d ; M 2d (C)) when 

£ ||p(*,y)|| 2 <oo. (3.4) 

(x,y)€Z d xZ d 

Note that (3.4) is equivalent to the Hilbert- Schmidt norm induced by the scalar product 
on l 2 {Z d x Z d ;M 2d (C)) 

(V, P) = Tr(7/*p) = Tr(r](x,y)*p(x,y)), (3.5) 

(x,y)& d x1 d 
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where we used the same symbol "Tr" for the trace in different spaces, which make £ 2 (Z d x 
Z rf ;M2d(C)) a Hilbert space. We also note that if p is non-negative, this implies for any 
7j d , (see [21] and Lemma 1.21 in [41]) 

p(x,y) = p(y,x)*, p(x,x)>0, and \\p(x, y)\\ < \\p(x, x)\\^ 2 \\p(y, y)\\ 1/2 . (3.6) 

Thus VM and the set of non-negative trace-class operators belong to l 2 (Z d x Z d ; M 2 d(C)). 

If po denotes the initial density matrix, its evolution at time n under U (n, 0) defined by 
(2.5) is given by 

p n = U(n,0)p U*(n,0). (3.7) 

The kernel of p n reads 

p n (x,y)= Yl J x k - k (n)po(x-k,y-k')J y k r k '*(n), (3.8) 

(k,k')e2. d xZ d 

and the expectation of the lattice observable F = I (g) / is denoted by 

(F) po (n) = Tr(p n (I®f))= £ Tr(p n (x,x))f(x), (3.9) 

if it exists. Again, we can express (F) po (n) as the expectation of a random variable on the 
lattice Z d : 

(F) po (n) = E w{n) (f(X n )), with Prob(X n = k) = W k (n) = Tr(p n (k,k)). (3.10) 

In case the evolution is random, the distribution W^^n) is random and the density matrix 
p% is random as well. We consider thus the expectation with respect to the randomness, 
noted E of the quantum mechanical expectation of the lattice observable, i.e. 

= E Wu{n) (f(X n )) = E w(n) (f(X n )), (3.11) 

where the distribution w(n) on Z rf is given by 

Prob(X n = k) = w k (n) = E(WHn)) = ETr(p^(A:, k)). (3.12) 

The corresponding characteristic function is defined by 

K°(y)=^ w (n)(e iyX ") = £ jy*Tr(E((%(x,x))). (3.13) 

The following assumption gives the required regularity properties on the lattice observ- 
able F = I <g) / and the initial density matrix po to legitimate the manipulations that 
follow. 

Assumption R: 

a) The lattice observable is such that, for any jjl < oo, 3C P < oo such that 

\f(x + y)\<C p \f(x)\, V(x,y) £ Z d x Z d with||y|| < p. (3.14) 



8 



b) The kernel po{x, y) is such that 

\\P0&V)\\ <°° (3-15) 

(x,y)<=Z d xZ d 

\f(x)\\\po(x,x)\\ <oo. (3.16) 

xez d 

Lemma 2.11 of [21] applies here as well, with the same proof, to ensure that for any n G N, 
the kernel p n (x,y) satisfies Assumption R if the kernel po does. For more discussion about 
properties of the density matrices we refer to [21]. 

3.1 Feynmann-Kac-Pillet formula 

We denote by I 2 (CI; M 2d (C)) the finite dimensional Hilbert space of M2d(C)-valued functions 
defined on Cl with scalar product defined by 

(<p,^) =Y,P(vm<P*(riMri)), (3-17) 

where the measure p on Q is the stationary initial distribution. We denote by |t)(t'| G 
l 2 (fl; M 2d (C)) the constant map which assigns |t)(t'| to any rj G Q and stress that the r, r' 
element of a matrix p G M 2d (C), r, t' G I±, can be expressed as 

(p) T y = Tr(\r')(r\p) = Tr((\r)(T'\y p). (3.18) 

Consider now the extended Hilbert space l 2 (1 d x 1 d ; l 2 (C>; M 2d (C))) ~ l 2 (1 d x Z d x 
$7; M 2rf (C)). Any p G / 2 (Z d x Z d x $7; M 2d (C)) can be expressed as 

p= (p(x,y;r]))^ y . rt)eZ d 

xZ d xfli 

where p(x, y; rj) G M 2d (C) (3.19) 

satisfies 

Y p(ri) TT (p( x ,y,v)*p(x,y,v)) < oo. (3.20) 

veil 

(x,y)GZ d xZ d 

The following is a version of Feynman-Kac-Pillet formula in the current setting. Let 
Po G l 2 (jL d x Z rf ;M2rf(C)) denote the initial density matrix, its evolution at time n under 
the random evolution operator U(n,0) defined by (2.5) and (2.6) is given by 

p n = U(n,0)p o U*(n,Q). (3.21) 

Since l 2 (Z d x Z d ;M 2d (Cj) <-> l 2 (Z d x 1 d x Q;M 2d (C)), we can consider p an element of 
l 2 (Z d ®Z d x ft; M 2d (C)), keeping the same notation. With the notation S x = \x) we have 

Proposition 3.2 Let K = l 2 (1 d x 1 d x Cl;M 2d (C)) and assume C holds. Then, if p G K, 
we have for any n G N, and any r, t' £ I±, 

E(p n (x,y)) T y = (5 x ®5 y ®\T)(T'\,M n p ) K , (3.22) 



9 



where the single step operator M : K, — ^ K, is given by 

(Mp)(x,y;r l ) = Y, Q(v,C)PA^r(r)V)p^-r(T),y-r(T'),C)^y-r(r')VTPr', (3.23) 

T,r'e/± 

(en 

where p G l 2 (Z d x Z d x Q; M 2d (C)) and Q(r?, C) = Prob(w(0) = r/|w(l) = Q. 
Remarks 3.3 i) Using that the initial distribution is stationary, it is easy to see that 

Q(C,r?) = ^|p(r ? ,C). (3.24) 
ii) In view of (3.12), the averaged distribution w(n) reads 

Wx{n) = Y E (pn(x, x)) TiT = (V x , M n p ) where V x = 5 x ®5 x ®Id. (3.25) 

rel± 

Hi) The adjoint of M, M* , acts as follows 

(M*p)(x,y; V ) = ]T Pfa, Q{o x Q*P T p{x + r(r), y + r(r'), QPA°vO- (3-26) 

r,r's/ ± 

iv) If {p(x,y,r])} X! y eZ d is self-adjoint, the same is true for {(Mp)(x, y; f])\ x y ^ I d . Such 
initial conditions p yield real valued quantities w x (n) = (vD^, M ra po)- 

Proof: First note that 

(S x ®5 y ® \t)(t'\ , M n Po ) !C = J>(C)((M>o)(x, y; ()) r/ . (3.27) 

(en 

Let £j = 53s=i r ( T s) an d ^ = Ss=i r ( r s)- Using the definition of M, we see that 

(M>o)(^,y;C)= Q((,Vn)Q(vn,rin-i)---®(m,vi)x 

{t(,-,<}6/„ 

P Tn (<Tx-tnO P Tn-l (^X-tn-lVn) " " " P T! (<Tx-tim)M X ~ *1) 2/ ~ *l) X 

(o- y _ t[ ri 2 yp Ti ■ ■ ■ K-c^)*^-!^-*^)*^- 

Since the initial distribution p is stationary, a straightforward computation shows that 

P(C)Q(C, ^)Q(%, ^-1) • • • Qfe, r/i) = p(r/i)P(?7i, r? 2 ) ■ ■ ■ P(ffr-i, »fc)P(%, 0- (3-28) 
Therefore, 

(<5^<5j / ^|r)(T , |,M n / 9 >^= »&) ■ ■ ■ P(??n-1, ??n)P(??n, C) x 

{•""!,■•■ ,T n }ei n 
M,-,rA}e/n 

{ll,12,' 1 1n,(}6!!"+ 1 

{APrn^x-tnQPrn-A^x-tn-irin) ■ ■ ■ Pn (<T x -ti V2) Po {X ~h,y- t[) X 

(o-y-t^mTPrl ■ ■ ■ ^y-t'^rinTPr^y-t'STPrk A ( 3 ' 29 ) 
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On the other hand, 

E(pn(z, y)) = ^2 Prob (fx-M^ii) = rn,a y _ t <u(i) = 77- for all z G {1, ■ • • , n}) 

{ti,--- ,T„}eIn 
{•>-{,■■■ ,^}e/„ 

{'71,'72.'.'7n}eri n 

X PrnVnPrn-iVn-1 • • • P^VlPO {x ~ h, y - t^) rf* ■ ■ ■ P T ' n _^in P r' n ■ (3-30) 

However, it is easy to see that 

Prob (a x _ ti u;(i) = r]i,a y _ t >u(i) = ??• for all i G {1, • • • ,n}) 

= f Prob ( W W = °x-uVi for all i € {1, • ■ ■ , n}) if <r~i u fH = %l t ^i ^ ^ 

otherwise. 

Now by letting aj = a~^ t .rji, and using that a; is a Markov chain on £1 we get 

E(p n (x,y))= J2 p(ao)P(o ,ai)P(ai,ta 2 )---P(a n _i,a n )x (3.32) 

{r(, 

{a ,a 1 r,an}en n + 1 

PrniVx-tna^Prn-! " " " P n (c^-ti «1 )P0 (x ~ *1 , 2/ ~ X 

K_ t /Q;i)*P T / • • • ^^K-t^i)*^;- 
Comparing (3.29) and (3.32) completes the proof. H 

3.2 Spectral Analysis of M 

Using Feynman-Kac-Pillet formula, studying the time evolution of our systems relies on 
the spectral analysis of the "single-step" operator M defined in (3.23). In order to do that 
we first take a closer look at the underlying symmetries of the systems. The operator M 
commutes with a group Q of unitary operators generated by translations: 

1. Simultaneous translation of position and disorder by an arbitrary element £ of Z d : 

<% p(x, y; u) =p(x-£,y- £; a^u), 

2. For r\ G T C II 1 such that a v = Id , M commutes with translation of the first position 
coordinate by rj: 

S^p(x,y;u) = p(x - n,y;uj). 

Note that S^S^ = S^, so the group of symmetries Q is isomorphic to Z rf x 7L d . We 
have chosen to use translation of the first position in the definition of S^; however, since 
a v = Id, we have S T] S^p(x, y; uj) = p(x, y — rj; oj). 

Remark 3.4 For any r\ G T = {£ G 7L d : = Id}, and any x G 7j d , we have a x + v = o~ x . 
Moreover, for any x G Z d , there exists a unique xq G Br such that a x = a Xo . 
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In order to take these symmetries into account in the spectral analysis of M, we define 
a generalized Fourier transform similar to [18]. We shall use the following notations 

L 2 (X;M 2d (C)) = {f:X^ M 2d (C) : ||/|| 2 = Tr(f /) = / dm(x)Tr(f* (x)f(x)) < oc}, 

Jx 

where m is a locally finite positive measure on X. 

Also, we introduce T* = {p* G M d \ p*j G Z, V7 G T}. If {7j}^ =1 is a basis of T, 
let {p)} d j= i be the basis of T* defined by p*ji = 5ij. We have that Z d C T*. We note 
T d = {p = Y?j=iPjP*j I Pj G [0,2vr),j = l,...,d}. With the linear map P : M d -»• M ' 
defined by its action on the vectors of the canonical basis as Pej = p*j, j = 1, . . . , d, we 
have Tp = PT d where T d denotes the d-dimensional torus. In particular, for any / defined 
on 

/ f(p)dp= f f(Pt)\P\dt, and Vol(T r d ) = (2vr) d |P|, (3.33) 

where | • | denotes the Jacobian determinant here. We denote the normalized measure on 
by dp = dp/{(27r) d \P\). 

We are ready to introduce the map F : l 2 (Z d x Z d x Q; M 2d (C)) -> L 2 (B r x T d x Tp x 
Q;M 2d (C)), defined by 

(W)(x, k,p; C) := *(x, k,p; C) = £ e wi-v)-ik-^ {x _ £ _ ^ _^ {3M) 

ver 

Since we can add to x any vector of T without changing the RHS, actually depends 
on xq G By defined according to remark 3.4. One checks that this generalized Fourier 
transform is a unitary operator with inverse 

(T- 1 X )(x,y;C)= I e- i y k e-^ x -^ x (x-y,k,p;a y C)dkdp, (3.35) 

JT d xT d 

where dkdp is the normalized measure on T d x Tp. 

Remarks 3.5 i) If {F^){x, k,p; Q = ^(xq, k,p; (), it satisfies for any p* G T* , any rj G T 

and any k* G Z d 

%o + v, k + 2irk*,p + 2tt P *; Q = e i27rp ^ $(x , k,p; Q. (3.36) 

ii) The operator {ip(x,y, ()} x ,y& d ^ s self-adjoint, i.e. ip(x,y,()* = ip(y, x ,() if and only if 
$(x ,k,p;() = $((-x) ,-k,p- k;a X0 Q*. 

Because of the symmetries of M, its expression TMT^ 1 in Fourier space admits a fiber 
decomposition of the form 

re 

7MT- 1 = / M(k,p)dkdp, (3.37) 

JT d xT d 
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where M(k, p) is an operator on I 2 (£?r x 0; M 2 d(C)) which becomes a multiplication operator 
in the variables (k,p) G T rf x Tp which we compute. The following expression holds for the 
(k,p) dependent "single-step" operator M(k,p) on L 2 (B r xT d xTfx Q;M 2d (C)): 

(M{k,p)^){x,k,p; V ) = Q(»7,C)e ifer(T V p(r(T) - r(T ' )) x (3.38) 

ceil 

Pr((r x -r(r)V)^ ~ t{t) + t{t) , k, p, (J- r {r')C) ((J-r(r')V)* Pr' ■ 

Remark 3.6 The action of the adjoint of M(k,p), denoted by M(k,p)*, reads 

(M(k,p)**)(x,k,p;ri) = P(r/,C)e- ifcr{r ' ) e- ip(r(T) - r(r ' )) x (3.39) 

r,r's/ ± 

cen 

(o- x t)*P T V (x + r(r) - r(T'),k,p,a r{Tl) () P T ,Q. 

Let us now consider the operator M(k,p) for (k,p) G T d x T r fixed as an operator 
on l 2 (Br x J);M2d(C)). As £?r and Q are finite, M(k,p) can be represented by a square 
matrix of dimension 4d|Br||0| depending parametrically on (k,p). Moreover, the map 
(k,p) ' ^ M(k,p) is analytic on C rf x C d . We denote the norm on l 2 (B r x fi;M 2 d(C)) by 

Proposition 3.7 Lei Spr denote the spectral radius, then for all (k,p) G T d x Tp 

Spr(M(fc,p)) < ||M(A;,p)||, a < 1, (3.40) 
On i/ie oi/ier /iond /or = and all p G Tp 

Spr(iQ r (0,p)) = ||M(0,p)|| /2 = l, (3.41) 

Remark 3.8 follows that Spr (M) = ||M|| = 1, where M is viewed as an operator on 
L 2 (B r x T d x T r x 0; M 2d (C)). 

Proof: First note that M(k,p) can be written as M(k,p) = TiSQ, where 

{QV)(x,k,p;rj) = J]Q(»7,C)tf(3,fc,p;C) 
(en 

(SV)(x,k,p;ri) = (a x ri)^(x,k,p;rj)(a vT (3.42) 

We fix (fe,p) and consider these operators on Z 2 (-Br x M 2d (C)) . Now Q = Id<8>Q where 
Q : l 2 (n;C) ->• Z 2 (ft;C) given by (5/(77) = £ C6 n Q(r/, C)/(C) and Id means the identity on 
l 2 (Br; M 2d (C)). An easy calculation using Jensen's inequality shows that for all / G / 2 (fi; C) 

2 

iiq/ii 2 = I>(r?) 



Ceo 



<E^)E^'0l/(C)| 2 (3.43) 
= E^)I/(C)| 2 = II/I| 2 . (3.44) 
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Therefore, we have that \\Q\\p < 1. 

On the other hand, for all * G l 2 (B r x 0; M 2d (C)) 



||s*||f a = 2p(0^[C**(^;^)(^C)*KC)*(^»/)C1 

i6B r 
C6S2 

= ^p(0^[**(x-,riMz;v)] = \\np- (3-45) 

Where we used the cyclicity of the trace and that elements of Q are unitary matrices. 
Finally to see that for any (k,p) G T d x Hp, ||S|| = 1, we notice that 

TV[(E*)*(x;C)(E*)(x;C)]= (3.46) 
E (a|**(x - r(r) + r(a); cr„ r(a) C)r)(r|^(x - r(r) + r(a); cr_ r(a) C)a) 

Now for fixed a, r let y = x — r(r) + r(a) and 77 = a_ r ^(. Using that a x are measure 
preserving transformations on S7, we have 

\\zml= E E^)H**(!'i'/M(^(i'i'/)«) ( 3 - 47 ) 

= J] p(r/)Tr r?)*^; r?)] = ||*|| 2 . (3.48) 

j/esp 

Putting the estimates on the norms of Q, S and S together we get the required bound on 
the norm of M(k,p) for all (fc,p) £T d xTf. 

Now consider £) = #o ® Id, where Id 6 / 2 (S7; M 2 ^(C)) takes the constant value Id. 

We compute 



(M(k,p)^i)(x; rj) = ^r{r') ^p{r{r)-r{r')) ^ q 



T, T 'ei± 



Pr^x-r(T)V){^-r(r')VT Pr'M X ~ K T ) + r(r')) 

= £ e ^( T ') e i f( r W- r ( T '))p T 5 TiT , ( 5o(x-r(T) + r(r / )) 

T,r'el± 

= E e ifcr(r) P T (5o(a;). (3.49) 



r£/± 



where we use X<efi Q( r ? ) C) = 1- From this it is clear that ^1 is an eigenvector of M(0,p) 
with an eigenvalue 1 for all p G Tp. Therefore we have that 

Spr(M(0,p)) = ||M(0,p)||,2 = l. (3.50) 
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Remark 3.9 A similar computations shows that 

M(0,p)*$i = $1. (3.51) 
If^i is considered a vector of L 2 (Br xT d xTf x Jl;M 2 d(C)) it corresponds to 

^i{x,y; V ) = ^- 1 $ 1 (x,y;r ? ) = <5 (z) ® <5 (y) ® Id ~ ld®|0)(0|. (3.52) 
^4/so ; to^/i t/ie definition (3.25) 

fc,p, r?) = e ife %(:r) Id = e ifc£ §i(x). (3.53) 

At this point we note that the characteristic function &n°(y) of the distribution w(n) 
satisfies, see (3.25) and (3.53) 

K°(V) = E e<WX <**> M >°> = E e ^<*- ^ ■) n )5b> /c (3-54) 

x<=Z d T 

In other words, slightly abusing notations, 

= ((M(y, ■)*)% Id, Po(y)} L2(BrxT , xfW2d(c)) 

= j jd (0*° M ' ^^)) n ^(y^)>P(B r xn;M 2d (C))^ 

Vpfa) / Tt(M(y,p) n i5b)(0,j/,p^)dp (3.55) 



where 



p (x, fc.p, C)=E e^-^-^/OoCx - i - V , -0 = Po(x, k,p) (3.56) 
>7er 

is independent of £• 
Remark 3.10 // 

Po = \ ( Po)( ( Po\ ® |0)(0| ~ Sq 5q ® | c,po) (v^o 1 5 where ipo G C d is normalized, (3.57) 

i/ten 

PoOej ? ?) = <5o(^) ® |v 5 o)(v 3 o| := Ro(x) is independent of (k,p,r]) (3.58) 

and 

$r?(y) = E^)/ TT(M(y,p) n \w)(<p \)(0,y,p, V )dp. (3.59) 
„en r 
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Hence, in the diffusive scaling, we need to control the large n behavior of the vectors 
M* (y/^/n,p) n 5 (g>ld and Mv/V™) in L 2 (B r xTfxU; M 2d (C)). This can be done, following 
the arguments of [21], under some spectral hypothesis. We shall discuss the validity of this 
hypothesis for specific cases later on and proceed by showing that it is sufficient to prove 
the diffusive character of the (averaged) dynamics, arguing as in [21]. We shall refrain from 
spelling out all details, referring the reader to the above mentioned paper. 

We work under the following spectral hypothesis on the matrix M(0,p) on l 2 (B^ x 
f2; M 2f j(C)). Let D(z, r) be the open disc of radius r > centered at z G C. 

Assumption S: For all p G T d , we have 

a(M(0,p)) n dD(0, 1) = {1} and this eigenvalue is simple. (3.60) 

Remark 3.11 Actually, because of (3.55), it is enough to assume that M(0,p))\x satisfies 
assumption S, where X is the M(k,p)* -cyclic subspace generated by 5q <g> Id, for all (k,p) G 
T d x T r . 

By analytic perturbation theory, there exist 5 > 0, v{5) > and k(6) > such that for 
all (k,p) eBfx T,f, where B d K = {y G C d | ||y|| < k} and T d = {y = yi+iy2, \yi G Tf,y 2 G 
M. d with \\y2W < f} the following holds: 

(7(% lP ))nJ)(M) = {\ 1 (k, P )} 

a(M(k,p)) \ {Ai(fc,p)} C D(0, 1 - S), (3.61) 

and the eigenvalue X±(k,p) is simple. For such values of the parameters (k,p) we have the 
corresponding spectral decomposition 

M(k, p) = X 1 (k,p)P 1 (k, p) + Mp i (k, p) (3.62) 

where M Pl (k,p) = P 1 (k,p)M(k,p)P 1 (k,p) and P\(k,p) = I - Pi(k,p). 

The simple eigenvalue X±(k,p), the corresponding spectral projector P\{k,p) and the 
restriction Mp 1 (k,p) are all analytic on B d x T d and Spr (Mp 1 (k,p)) < 1 — 5. Moreover, 
for any p G Tp 

lim XAk,p) = 1, and lim Pi(fc,p) = |*o)(*o| = n > (3-63) 

k^O k^O 

where $ = ^i/||*i|| = ® Id, see (3.51). 

Taking into account the fact that w x (n) is real valued for any selfadjoint and trace class 
po, we have 

$«>(&) = $£°(-fc), for all k G T d . (3.64) 
This yields a symmetry of Ai: 

Lemma 3.12 For all k G B d , and all p G T r , the following identity holds 

Ai (k, p) = Xi{-k,p-k). (3.65) 
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Proof: It follows from (3.55) that for any k £ R d and any self adjoint trace class po 
K°(k) = ^ d < 5 o®Id, M(k,p) n Mk,p)) P{BrXn . M2dic)) dp 

The first step consists in showing the pointwise identity of the smooth scalar products in 
l 2 (B r x Q;M 2d {C)) for p = 5 (x) ® \<po){<po\ = R {x) : 

(6 ® Id, M(k,p) n R ) - (5 ® Id, M{-k,p - k) n R ) = 0. (3.67) 

Identity (3.66) holds for any self-adjoint po, thus in particular for po(xq, k,p) = b(k,p)Ro(x) 
where b belongs to the vector space of periodic functions satisfying 

b : T d x Tjf C, such that b(k,p) = b{-k,p- k), (3.68) 

see Remarks 3.5. Therefore, we get for any such b 

= / (5 (8) Id, M(k,p) n R )b(k,p) - (5 <g> Id, M(-k,p) n R )b(k,p + k)dp 

= I (<<5 <8» Id, M(k,p) n R ) - (5 ®Id, M(-k,p-k) n R )y(k,p)dp. (3.69) 

An example of smooth function b satisfying our requirements is b\{k,p) = fi(k)g±(2p — k) 
with 

gi : R d -> R, 2tiT* - periodic, and /1 : R d -»• R, 2vrZ d - periodic and even. (3.70) 

Note that Z d C T* ensures Z periodicity of 61 in and that b\(k,p + 2-7T7*/2) = b\(k,p), 
for all and 7* 6 T*. Another slightly more complicated choice is constructed with 

g 2 : R d -> R, 27rr* - anti-periodic i.e. # 2 (x + 27r 7 *) = -pa (a:), Vi = 1, . . . , d, (3.71) 

and {7* } d =l is the basis of V* . Then f 2 : M d ->■ R is defined as follows: for j = l,...,d, 
e j = 52i=i m i(jhh where rrii(j) G Z, and 

/ 2 (x + 27re j ) = (-l) E '= imi(j ' ) /2(a:), VxeR d , Vj = l,...,d. (3.72) 

That is f 2 is 27r-periodic or 2-7r-anti-periodic in the direction ej, depending on the compo- 
nents of the corresponding vector 7*. Then, by construction, b 2 (k,p) = f 2 {k)g 2 (2p — k) 
satisfies our requirements and, moreover b 2 (k,p + 2ti^*-/2) = —b 2 (k,p). 

Now, assume (3.67) doesn't hold at some po £ Tp. By a suitable choice of g\ and g 2 as 
above, we can construct a smooth b(k,p) = b\(k,p) + b 2 (k,p) that is non-zero in a small 
neighborhood of po only so that (3.69) fails, which yields a contradiction. 
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Then one exploits the spectral decomposition (3.62) and (3.63) with (8q ® Id, Rq) = 1 
to deduce from the above that if \\k\\ is small enough 

/ . — . n \ 1/ Tt 

\i(k,p) = lim (<5 <8> Id, M(k,p) R )\ 

1/n 



= lim (<5 ®Id, M(-k,p- k) R )) = X 1 (-k,p-k). (3.73) 

The result extends to complex k by analyticity of Ai(-,p) in M 

We now compute a second order expansion of X\(k,p) = Tv(Pi(k,p)M(k,p)) around 
k = 0, using the decomposition (3.42) 

M{k,p) =T,(k)SQ, (3.74) 

where only the unitary map £ depends on k (and p), as stressed in the notation. We expand 
£(/c) as 

E(k) = E(0) + Ei (A;) + E 2 (fc) + O p (||fc|| 3 ), (3.75) 

where Ej(/c) is of order j = 1,2 in k and the remainder is O p (||A;|| 3 ), uniformly in p G T„ ■ 
Explicitly, 

(£i(fc) + £ 2 (fc))tf(x,fc,p;C)= (3.76) 
^ (^r(r') - (A:r(r / )) 2 /2)e ip{r(T) " r(r ' )) P T ^(x - r(r) + r(r'), k,p; a _ r{jl) QP T , . 

Then, in terms of the unperturbed reduced resolvent S p (z) defined for any p G Tjf and 2 in 
a neighborhood of 1, by 

(M(0,p) - zr 1 = + S p (z) (3.77) 

we have see [23], p. 79, 

Ai(fc,p) = l + Tr(Ei(A;)5Qn) (3.78) 
+ Tr(E 2 (fc)5gn - Ei(fc)5Q,Sp(l)Si(fc)5Qn) + O p (||fc|| 3 ). 

Explicit computations making use of SQ^ = *o, SQS P (Y) = E(0) _1 (I - II + S p (l)) and 

(E(0)- 1 $)(x,r / )= ^ e^M-^'»P r cI>((x + r(r)-r(r / ), C 7 r ( T0 )P T ' (3.79) 

(T,r')e4 

yield 

Lemma 3.13 For all p G 7^ and G ,B(0, /c) ; i/iere exists a symmetric matrix D(p) G 
Al (A;,p) = l + ^ J>r(r) + P (IW| 3 ) 

r6/± 

= 1 + i E fcr ^ " |<fc|»Cp)*> + Opdlfef). (3-80) 
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The map p i— >■ D(p) is analytic on T d ; when p G Tp, £ M2d(M) is non-negative and 

n (p)iJ = dkfdkj ^'P^' G {1,2, ... Moreover, O p (||/c|| 3 ) is uniform in p G 7^. 

Proof: Existence and analyticity in p of D(p) follow from analyticity of Ai in y and 
analyticity of S p (l) in p, see (3.77). Since U)(p)ij = A(0,p), the matrix is symmetric. 

For p G Tp, Lemma 3.12 implies that A(0,p) is real valued, hence the matrix elements H>(p) 
for p G Tp are real as well. Finally, (3.40) implies that (k\B(p)k) > for all k G T d . ■ 

As a consequence of the spectral analysis above, it follows exactly as in [21], that 

Proposition 3.14 Under assumption S, uniformly in p G T d , in k in compact sets of C d 
and in t in compact sets ofM.+ , 

lim M{k/n,p)W = e ityT U, (3.81) 

lim M{k/^p)^e~ i[tn]Tyl ^ = e-^vWP^U. (3.82) 



4 Diffusion Properties 

These technical results lead to the main results of this section which are the existence of 
a diffusion matrix and central limit type behaviors in the diffusive scaling, as in [21], with 
the same proofs, that we don't need to repeat. 

Let A/"(0, X) denote the centered normal law in R d with positive definite covariance 
matrix £ and let us write X w ~ A/"(0, X) a random vector X w G M. d with distribution 
Af(0, £). The superscript co can be thought of as a vector in M. d such that for any Borel set 
A C R d 

F(X" G A) = 1 , I e-3< w l E_lw >dw. (4.1) 

V 7 (2vr) d / 2 7det(S) 7a 

The corresponding characteristic function is (y) = E(e* J/X ") = e~2^l s y). 

The first result concerning the asymptotics of the random variable X n reads as follows 
for an initial density matrix of the form po = | c^o) ( I <8> |0)(0| : 

Theorem 4.1 Under Assumption C and S, uniformly in y in compact sets of C d and in t 
in compact sets ofW + , 

lim ^l (y/n) = e it ^ (4.2) 
lim e- iltn] %*'« , ] {y/y/n) = f e'^ ^ dp, (4.3) 

where the right hand side admits an analytic continuation in (t,y) G C x C d . In particular, 
for any (i,j) G {1,2, . . . ,d} 2 , 

lim = n (4.4) 

lim ((X-nrUX-nr) j U (n) = f ^. {p)d ~_ (45) 
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IfU)(p) = ID > is independent of p 6 Tp, then, for any initial vector = ¥0 ® |0), we 
Ziaue as n -> 00, wi/t convergence in law, 



X n -nr v 



X" ~JV(0,D). (4.6) 



Remark 4.2 JFe u>iZZ ca// diffusion matrices &oi/i D(p) and D = J Td D(p) dp. 

Remark 4.3 VFe prove freZow f/iaf a central limit theorem for X n may hold in cases where 
D depends on p, see Theorem 6.4- 

For initial conditions corresponding to a density matrix po, we have 

Corollary 4.4 Under Assumptions C , S and R for the observable X 2 , we have for any 
t>0,y£C d , 



lim &L(y/n) = e*^ , (4.7) 



lim e V™ 



*&(^ = / e^«*^i|npo(-,o lP ,)} i2(BrxT , xQ;MM(c)) ^ 

•/lip 



e -ife|B(p)y> Tr(p o )(0,0,p)dp (4.8) 



where, see (3.56), 

p o (0,0,p) = ^e-^(e-(,0. (4.9) 



/or any (i,j) G {1,2, ... ,d} 2 , 



lim^)=F 1 , (4.10) 

n— >oo 7j 

lim «*-"rMX-nr) 3 -Un) = /" ^ ^(pbXO, 0>p) ^ (4 . n) 

From Corollary 4.4, and Theorem 4.1, we gather that the characteristic function of the 
centered variable X n — nf in the diffusive scaling T = nt, Y = y/y/n, where n — > 00, 
converges to 

e -£<.|D-i(p).> \ 



Tf ^((i^v™J fa)Tr(?0)(0 ' 0,p) * ( " 2) 

where the function under the Fourier transform symbol J 7 is a solution to the diffusion 
equation 

i,j=l % 1 

As explained in [22], [18], it follows that the position space density Wk([nt])5(y/nx — k) con- 
verges in the sense of distributions to a superposition of solutions to the diffusion equations 
(4.13) as n ->■ 00. 
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5 Moderate Deviations 



It is shown in [21] that the spectral properties of the matrix M(k,p) proven in Section 
3.2 allow us to obtain further results on the behavior with n of the distribution of the 
random variable X n defined by (3.12) with localized initial condition po = l^oKvol <8> |0)(0|, 
corresponding to the vector Rq G l 2 (Br x Q; M2d(C)) , see (3.58). This section is devoted 
to establishing some moderate deviations results on the centered random variable X n — nf. 
Again, since all proofs are identical to those given in [21], we merely state the results. 

Moderate deviations results depend on asymptotic behaviors in different regimes of the 
logarithmic generating function of X n — nf defined for y G M. d by 

A n (y) = ln(E w(n) (e^--))) G (-00,00]. (5.1) 

This function A n is convex and A n (0) = 0. 

Let {a n }neN be a positive valued sequence such that 

lim a n = 00, and lim a n /n = 0. (5-2) 
n— >oo n— >oo 

Define Y n = (X n - nf) / y/na n and, for any y G R d , let A n (y) = ln(E w ( n )(e yYn )) be the 
logarithmic generating function of Y n . 

Proposition 5.1 Assume C and S and further suppose D(p) > for all p G T d . Let 

y G R d \ {0} and assume the real analytic map T d B p \— > (y\0(p)y) G 1R+ is either constant 
or admits a finite set {pj{y)}j=i, -,j of non- degenerate maximum points in T d . Then, for 
any y G R d , 

lim — A n (a n y) = J(y|B( Pl (y))y) (5.3) 

which is a smooth convex function of y. 

Let us introduce a few more definitions and notations. A rate function I is a lower 
semicontinuous map from M. d to [0, 00] s.t. for all a > 0, the level sets {x \ I(x) < a} are 
closed. When the level sets are compact, the rate function / is called good. For any set 
r C M. d , r° denotes the interior of T, while T denotes its closure. 

As a direct consequence of Gartner-Ellis Theorem, see [16] Section 2.3, we get 

Theorem 5.2 Define A*(x) = sup^d ((y\x) - \{y\H(p 1 {y))y)) , for all x G R d . Then, A* 
is a good rate function and, for any positive valued sequence {a ra } rag N satisfying (5.2) and 
all Borel set V C R d 

- inf A*(x) < liminf — ln(P(Of n - nf) G Jn~b^T)) 

< limsup — \n(¥((X n - nf) G ^/na^T)) < - inf A*(x). (5.4) 
Remark 5.3 As a particular case, when D(p) = ID > is constant, we get 

A*(x) = )^(x\B- l x). (5.5) 
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Remark 5.4 Specializing the sequence {a n } ne N to a power law, i.e. taking a n = n a , we 
can express the content of Theorem 5.2 in an informal way as follows. For < a < 1, 

¥((X n - nr) G n( a+1 )/ 2 T) ~ e^^^^'W. (5.6) 

For a close to zero, we get results compatible with the central limit theorem and for a close 
to one, we get results compatible with those obtained from a large deviation principle. 



6 Large Deviations 

In this section, we push further the analysis of the large n behavior distribution of the 
random variable X n (defined by (3.12) with localized initial condition po = |v?o)(<A)| ® 
|0)(0|) by proving large deviations estimates and a central limit theorem under stronger 
assumptions on the spectral properties of the matrix M(k,p). 

We change scales and define for n G N* and y G R d a rescaled random variable and the 
corresponding convex logarithmic generating function 

Y n = Xn ~ nr and A n (y) = lnE w(n) (e^) G (-oo, oo]. (6.1) 

Because of the new scale n, the existence of lim^oo = lim^oo lnEw< - n ^- e 1 _ yf i s 

not granted for all y, by contrast with the previous section. However, because ||1^|| < cq, 
for some cq < oo, we have for any y G C d , and a fortiori for any y G M d , 

< bile (0-2) 

n 

Moreover, as the next Proposition states, the limit exists for ||y|| small enough, under more 
global, yet reasonable, hypotheses: 

Proposition 6.1 Let y G M. d n B(0,k) be fixed and assume the function Tf B p t— > 
|Ai(— iy,p)\ is either constant or admits a finite set of non- degenerate global maxima {pj(y)}j= 
in Tp. Further assume V p \i(—iy,pj(y)) = 0, for all j = 1, . . . , N. Then, for k > small 
enough, 

\(ny) 

lim = -yf + ln(\ 1 (-iy,pi(y))) (6.3) 

n-¥oo n 

is a smooth real valued convex function of y G B(0, k) n M. d . 

Remarks 6.2 i) In case Xi(—iy,p) = \\(—iy,0) is independent of p G T d , the right hand 
side of (6.3) equals — yf + ln(Ai (—iy, 0)). 

ii) The assumption V p Ai(— iy,Pj(y)) = may be too strong to deal with certain cases. 
However, if it does not hold, in which case V p Ai(— iy,Pj(y)) G iM. d , the asymptotics of the 
integral that yields A(ny)/n is out of reach of a steepest descent method without further 
information on the behavior of X±(—iy,p) for p away from Tp. 
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The proof is a straightforward alteration of that of Proposition 5.1, based on Laplace's 
method to evaluate the asymptotics of the integral 



M 2d (C)) 



dp (6.4) 



Xe nin(Ai( -^ )} <^ 



where 7 > and the prefactor is non-zero, due to the smallness of ||y||. For complete- 
ness, we briefly recall the argument in case there is only one maximum at p\ 6 Tp. 
Dropping the variable y in the notation and writing ln(Ai(p)) = a(p) + ib(p), V(p) = 
(tyi\P(— WiP)Ro)i2( BrX Q M 2d (C)) we nave i n a neighborhood of p\ G Tp determined by 
Va(pi) = and D 2 a(p 1 ) < 

e nln(Ai(p))p(p) = e ™MAi(pi))p( pi ) x ( 6 .5) 
xe in(Vb(p l )(p^p l )) e n((p^p l )\(D 2 a(p 1 )+iD 2 b( Pl ))( P - Pl ))/2 e nO(\\p-p 1 \\'- i )^ + - pi||)). 

Making use of D 2 a{p\) < 0, we can restrict the integration range in (6.4) to B(pi,fi(n)) C 
M d , with l/^/n <C /Lt(n) <C 1/n 1 / 3 , at the cost of an error of order e~ nfM ^ c , for some c > 0, 
so that we are lead to 

l e m(Vb( Pl )b> e n( P |(D 2 a( P1 )+iD 2 6( Pl )) P >/2^ ^ + 0( n/i ( n )3) + 0(jl{n))). (6.6) 

JB(0,n(n)) 

When Vb(pi) 7^ 0, the analysis of the large n behavior of (6.4) and (6.6) requires global 
informations about the analytic properties of Ai for p far from the real set Tp, hence we 
require Vb(p\) = 0. Since Ai = 1 + 0(||y||) / 0, we have 

Vo(pi) = & KAi(pi)VKAi(pi) + 9Ai(pi)V9Ai( m ) = (6.7) 

V 6 (P.) = V« g (A l(p) ) = ^Ej>, (6.8) 

so that the hypothesis Vb{p\) 6 M rf implies VAi(pi) = 0. Now, at the cost of another error 
of order e~ n ^ n ^ c , (6.6) equals 



/ 



e n(p\(D a( P1 )+iD b( P1 ))p)/2 dp ^ + ( n/ ,( n )3) + 0( M (n))) + C ), (6.9) 

where a Gaussian integral yields 

/• e n( P \{D^) + iDH {p ,)) P )/2 d G where (2vr) rf / 2 
V n d / 2 ' ^det(- J D 2 a(pi) - iD%( Pl )) 

Altogether, we get 

ln(E U)(n) (e^)) 



n 



= ln(Ai(-iy,p)) (6.11) 
+ 1 In (^^(1 + 0(n^(n) 3 ) + 0( M (n)) + 0(e~ n ^ 2c )) , 
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which yields the result in the limit n — > oo. H 
We set for all y £ R d 

A(y) = limsup^^ G (-00,00), (6.12) 

which is convex, finite everywhere and bounded by co||y||. Moreover, for ||y|| < n, A(y) 
equals the right hand side of (6.3) and is thus smooth, and A(0) = A(0) = 0. Let us 
consider the Legendre transform of A 

T(x) = sup ((y\x) - A(y)) > 0, for all x G R d . (6.13) 

We are now in a position to state our large deviations results via Gartner-Ellis Theorem. 

Theorem 6.3 Assume the hypothese of Proposition 6.1 Let A and A be defined by (6.12) 
and (6.13). Further assume A is strictly convex in neighborhood of the origin. Then A is 
a good rate function and there exists rj > such that for any T 6 R d 

1 — * 

limsup-ln(P((A n - rif) G nT)) < — inf A (x) (6.14) 

n xe r 

1 — * 

liminf - \n(¥((X n - rif) 6 nT)) > — inf A (x). (6.15) 

n^oo n xer°nB(o,r,) 

Proof: Exercise 2.3.25, p. 54 in [16], shows that since A is finite on M d then A* is a good 
rate function and that (6.14) holds. 

To show that (6.15) holds, we invoke Baldi's Theorem, Thm 4.5.20 in [16]. First, Exer- 
cise 4.1.10 of [16], point c), shows that the law of Y n is exponentially tight, as a consequence 
of A being a good rate function and (6.14) holding true. Then, by Exercise 2.3.25 still, 
if x = VA(y) = VA(y) for some y £ B(0, k), then x £ J 7 , where T is the set of exposed 
point for A with exposing hyperplane y. Let us recall that this means that for all z 7^ x, 
yx — A (x) > yz — A (z). Now, since A is strictly convex at the origin, its Hessian at zero 
is positive definite and VA(0) = 0. It thus follows from the implicit function theorem that 
for some rj > 0, the map y — > VA(y) is a bijection with range B(0,r]). Hence B(0,r]) is 
included in the set of exposed points for A . Also, the corresponding set of exposing hy- 
perplanes belongs to -B(0, k), where A coincides with A, which is finite everywhere. Hence, 
all hypotheses of Baldi's Theorem are met, so that (6.15) holds. ill 

Another direct consequence of Proposition 6.1 together with (6.2) is a central limit 
theorem for X n , as proven by Bryc, [13]. A vector valued version of Bryc's Theorem suited 
for our purpose can be found in [19]. 

Theorem 6.4 Under the assumptions of Proposition 6.1, we have, with convergence in 
law, 

^L_^!^AA( ,B) (6.16) 

where ®ij = aSr A (°) > 0. 
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Remark 6.5 The results of this section carry over to the cases considered in [21], see also 
Section 9. 



7 Example 

Let us consider here a fairly general situation in which the spectral hypotheses we need can 
be checked explicitly. 

We work in Z d and consider a model characterized by a representation of Z d , x i-> a x 
of measure preserving maps, a jump function r : I± — > 7j d such that 

r{T)-r{T') e T, Vr,r' G 7±, (7.1) 

a kernel P with identical lines 

Pfa,C) = P(C), Vt?gO, (7.2) 
and a set of unitary matrices {r]} v ^Q with trivial commutant 

WUn = {cI,c€C}. (7.3) 

This implies that the corresponding stationary distribution is p(Q = P(C)- We address the 
simplicity of the eigenvalue 1 of M(0,p)|x, see Remark 3.11. 

Proposition 7.1 Under assumptions (7.1), (7.2) and (7.3), M(k,p)\x is independent of p 
and M(0,p)\x admits 1 as a simple eigenvalue. 

Proof: The simplicity of the eigenvalue 1 of M(0,p)\x is equivalent to the simplicity of 
the eigenvalue 1 of M(0,p)*\x- 

We first observe that M(k,p)* leaves the subspace 

J = span{5 <8> A | A : tt ->■ M 2d (C) is constant} (7.4) 

invariant: 

(M(k,p)*5 ®A)(x,ri) (7.5) 

Ceil 

= e^5 (x) E p{Q(.v x Q*P T AP T ,e- ik ^'\ 

Ceo 

= <fo(x)X>(C)CMtf(*OC, 

fen 
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where U(k) = Yl T 'ei± ^ >T ' e lkr ^ T '^ ■ Hence we have X C J and M(k,p)*\j is independent of 
p G Tp. Thus we can consider M(0,p)*\j. 

Note that U(0) = I and that M(0,p)*\jS A = 5 (g> A is equivalent to M(A) = A 
where 

M(A):=Y J P(C)CA(, VieI 2i (C). (7.6) 
fen 

With the scalar product (A\B) = Tr(A*B) on M 2d (C) we have 

||X(^)|| 2 = ^ p(CM»/)(»/*^|CMC), (7.7) 
(C,»?)en 2 

where | (r/* J 4r/|C* J 4C)l < ll^4|| 2 > with equality if and only if if At] = e ie ^CAC, for some 
6 V<: G R. Hence ||7W(A)|| = ||A|| if and only if r]*Ar] = (*A(, for all 77, C- Thus, any 
invariant matrix under Ai satisfies 

M{A) = rj*Arj = A, V77GO. (7.8) 

Since the commutant of {r]} v en is assumed to be reduced to cl, c G C, we get the result. I 

8 Examples of diffusive random dynamics 

In this section we consider a specific example of measure dfi on U(2d), the set of coin 
matrices, for which we can prove convergence results on the random quantum dynamical 
system associated with (3.1) for large times. This example is a generalization of the example 
considered in [21] for site- independent coin matrices. While the following results hold for 
vector and density matrix initial conditions, we only consider here the vector case, for 
shortness. 

8.1 Permutation matrices 

We start by recalling a few deterministic facts. Let &2d be the set of permutations of the 
2d elements of I± = {±1, ±2, . . . , ±d}. For ir G &2d, define 

C(tt) = £ K(r))(r| G U{2d) so that C^tt) = S aMr) , (8.1) 

re/± 

and C(tt) is a permutation matrix associated with ir. Note the elementary properties: for 
any ir, a G 6 2 d, 

C(I) = I, C*(tt) = C t (tt) = (^(vr" 1 ), C(7r)C(<7) = C(na). (8.2) 

The matrices C(ir) allow for explicit computations of the relevant quantities introduced 
in Section 2. Given a sequence {Cj = C(iTj)}j = i t ... ;n of such matrices, a direct computation 
shows that with the definition Tj = 7Tj(rj_i), Jj^(n) takes the form 

4(n) = £ I^X^M, (8.3) 
rie/-]-. s.t. 



26 



and J°(n) = 0, if £™ =1 r(r,)/fc. 

Consequently, the non-zero probabilities Wk(n) on Z d read for any normalized internal 
state vector y?o- 

O) = Pl{n)M 2 = E IK 1 ^)!^)! 2 . (8.4) 

riG-Tj- s.t. 

Moreover, with ri = 7Ti(to) we get 

<f>0 = E aT ol r °) ^ K^r 1 ^!)!^)! 2 = E Ko| 2 <Wi(t )- (8-5) 

ro6/± r el± 

Hence °(n) = £ roe/± K | 2 <5 E « =1 rfc),* so that forf = I ®/ and = ® |0)<0| 

n 

(F)^ (n) = E («)/(*) = E KlV(I>fa)). ( 8 - 6 ) 

fcez d T () e/ ± j=i 

Remarks 8.1 In other words, given a set of n permutations, there is no more quantum 
randomness in the variable X n , except in the initial state. 

If one generalizes the set of matrices by adding phases to the matrix elements of the per- 
mutation matrices, it does not change the probability distribution {W£° (n)} feeZ d , see [21]. 

Therefore the characteristic functions take the form 

Corollary 8.2 With Tj = (iTjirj-i • • • 7Ti)(to), for j = 1, . . . ,n, 

*S°(v) = E e*£7=i r < T '>K| 2 (8.7) 

The dynamical information is contained in the sum S n = 5^j=i r ( T j) w hich appears in 
the phase. The next section is devoted to its study, in the random version of this model 
where the coin matrices are random variables with values in {C(n), ir G ©2d} distributed 
according to (3.1) 

8.2 Random Setup 

We consider that the permutation matrices are given by the process defined by (3.1) and 
we identify C(tt) and it: 

Assumption M: 

Let {o;(n)} ne N be a finite state space Markov chain on C &2d with transition matrix P 
and stationary initial distribution p and a representation a of 7L d of the form x ^ a x where 
for each x € a x : — > ft with C &2d, is a measure preserving bijection. We set 
C%{x) = a x {cj(n)) with C£(0) =u(n). 

We have that for every x 6 7j d , the set of random matrices/permutations {7r^(x))} n£ N = 
{a x (Lu(n))} n £fq, with u{n) Gflc ©2d, the Markov chain. 
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Given a set of random permutation matrices as above, we start at time zero on site G 
Z rf , with initial vector |to) <8> |0). The dynamics induced by the permutation matrices sends 
this state at time n > 1 to the state | YL™=i r ( T s)) ® \ T n), where Tj = a^y-i r ^ T ^(oj(j))Tj-\. 

Hence, in view of (8.7), we introduce the random variables S n (uJ) = J2j=i r ( r j(^)) e ^ 
and r(Tj(to)) where Tj(cJ) is defined for j = 1, . . . , n by 

ti(w) = <T (w(l))r , r,-(57) = a^-i^^Mj))^^), (8.8) 

for a given tq. Note that ^(cJ) = Tj((uj(j),uj(j — l),...,u;(l)). They have the following 
properties 

Lemma 8.3 Let (po = ^ T) ar | T o) be the initial vector, and assume M holds true. Let 
{Tj(aJ)}jgN be the I± valued process defined by (8.8). Then, with the notation 

Prob((r n (w),...,ri(w),r (w)) = (r n , . . . , n, r )) = T(r„, . . . , n, r ), n G N, (8.9) 
T(r n ,...,ri,To) = |a ro | 2 ^ p(vri)P(CT r{ri) (vri), vr 2 ) • • • P(o- r(Tn _ l) (vr n _i), vr n ) x 

7ri,...,7r n GSl 

x(r n |C(7r n )r n _i) • • • (n\C (m)^) . (8.10) 
Proof: We start with T(tq) = |a T0 | 2 ; according to the initial condition, and 
T(t!,t ) = |a T0 | 2 Prob(a;i s.t. 0o(wi)(t o ) = n) 

= l°ro| 2 *n,«»u(in)(7D)P( 7r i) = l a ro! 2 ( T i\ C (M*i)) r o)p(7Ti). (8.11) 
7rieSl 7rieS7 

Note that since oo is the identity, T(t\,tq) = ~E p ((ti\C (oj) ro))|a TO | 2 . Then 

T(t 2 ,ti,t ) = |a T J 2 Prob((a;i,a;2) s.t. <7o(wi)(r ) = n andov( n )(w 2 )(Ti) = r 2 ) (8.12) 

= l a r„| 2 ^2,^( T1 )( 7 r2)(T 1 )^r 1 , CT() (7r 1 )(r„)P(7ri)I !) (7ri,7r2) 
7ri,7T26n 

= Kol 2 ^ (^^(^(n)^)) n)(Ti|C(cr (vri)) T )p(7ri)P(7ri,7r 2 ), 

7ri,7r 2 eSl 

and, by induction 

T(r n , . . . ,ti,t ) = \a T0 | 2 ]T p(7ri)P(iri,ir 2 )---P(7r n _i,7r n )x (8.13) 

7ri,...,7r n eSl 

x(T n |C(<7 E n-i r(7i) (7r n ))r n _i) • • • (ti|C(o-o(7Ti)))to). 

Using the properties of a, the measure invariant representation of we get for any j > 1 
withTfj = a E ,-i r(rs) (^), 

2 P (^-l' 7r j)( 7 "jl C '( (7 ^-; r ( Ts )( 7r j)) r j-l) = X] P ( CT r(r J _ 1 )(^'-l),^)(^|C(^)T i „l), 

(8.14) 
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which ends the proof. 



The distribution of {Tj(uj)}j & ^ is neither stationary, nor Markovian, in general. But we 
can express it in a more convenient way as follows. 

Consider the space C 2d C' n ' with orthonormal basis denoted by {|r ® Tr)} T £i±,nen- 
Let N 6 M 2 d\n\Q& + ) be defined by its matrix elements 

{t'®w'\N t®tt) = {t'\C(tt') T )F(a r{T) (ir)y) = <V )7r , (T) P(<7 r(r) (7r),7r'), (8.15) 

and the vectors = ^ T ei± |t <g) 7r) and A(to) = Z^tt t l a T"ol 2 P( 7r )( T l^ Y ( 7r ) T o)l r 7r )- Then, 

7TSS2 

(8.10) reads 

T(r n , . . . ,ti, T ) = <*i|(|r n )(r n | ® I)iV(|T n _i)<T n _i| ® I) . . . AT(|ti)(ti| I)A(r )>. (8.16) 
Introducing also the matrices D(y) and iV(y) on C 2d (8) C'^' with y £ T rf by 

D( y ) = d(y) <g) I, where d(y) = ^ e iyr(T) |r)(r| and iV(y) = D(y)N (8.17) 

r£/± 

we can express the characteristic function <3?^ : T rf — > C of the random variable S n (uJ) = 
E^ir(r,-(5;))a8 

€(2/) = E e i «^=i r ^T(r n ,...,r 1 ,7 l) ) (8.18) 

Tn,T n -l,...,roG/± 

= (^KiVd/))"- 1 ^^)), where Bfo) = £ Afo). 

At this point, we can apply the same methods as above to describe the large n behavior 
of S n (uJ) by studying the asymptotic behavior of the suitably rescaled characteristic function 
&n(y) under appropriate spectral assumptions on TV. 

Note that N is a stochastic matrix and that P and p are invariant under a x so that we 
have 

iV T *i = *i, N X i = Xi, and || 7V|| = Spr (N) = 1 (8.19) 

with 

#i = E k® 71 ") and Xi = ^pW|t®7t). (8.20) 
Tren yren 
Also, .D(y) being unitary for y real, we have ||iV(y)|| < 1 for all y G T d . 
Assumption S: 

<j(N) n D(0, 1) = {1}, and this eigenvalue is simple. (8-21) 

Remarks 8.4 The corresponding spectral projector of N reads \xi)(&i\/(2d). 

Again, it is enough to assume that S holds for the restriction of N to the N(y)* -cyclic 

subspace generated by ^i. 
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The perturbative arguments given in Section 4 leading to Corollary 4.1 by means of 
Levy Theorem apply here. For y G C d in a neighborhood of the origin, let Ai(y) be the 
simple analytic eigenvalue of N(y) emanating from 1 at y = 0. Let v G M. d and the non 
negative matrix E G Md(R) defined by the expansion 

Ai(j/) = 1 + ij/tJ - ^(y|Sy) + 0(||y|| 3 ). (8.22) 

Explicit computations yield for any y € C d 



2d 



where 5(1) is the reduced resolvent of N at 1, 771 = 1 7r) and 771 = ^7r^( 7r )l 7r )- 

Proposition 8.5 Lei <£> = S T0 6/± a ^ol T o) anc ^ ^ S'n(w) = S"=i r ( r j(^)); w^/t Tj(uJ) 
defined by (8.8). Assume M and S and let E 6e defined by (8.23). Then, if E > , we 



/icwe as n — ^ 00 

© f (8 . 24) 

n 

S ^~ nl] ^ X<^Af(0,S). (8.25) 

V W 

As a consequence, for any sample of random coin matrices, we obtain the following long 
time asymptotics of the quantum mechanical random probability distribution W?°(n) of 
the variable X!f t , whose characteristic function is defined by (8.7). 

Theorem 8.6 Under the assumptions of Proposition 8.5, the following random variables 
converge in distribution as n — > 00: 

e -W^o (y/ ^ = J- \a T f (e iy ^ {Sn ^- n ^ — > e^, (8.26) 
r ei± 

where X w ~ A/"(0, E). Moreover, for any (i,j) G {1,2,..., d} 2 , as n — > 00, roe /iave m 
disiri&tiiion, 



n 



fi (8.27) 



((X-nr),(X-nr)^ o (n) ^ 
n ifc 
where B^ fc is distributed according to the law of X^X^, where X w ~A/"(0, E). 

Proof: Identical to that of Corollary 6.8 in [21]. 



(8.28) 
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8.3 Specific Case 

Let us close this section by providing an example that satisfies the assumption S. It is 
the case where the kernel P depends on the second index only, i.e., when the permutations 
{^COIjeN are i.i.d. and distributed according to p. 

Proposition 8.7 Assume M with a kernel P satisfying P(tt',tt) = p(ir). Let P be the bi- 

stochastic matrix acting on C 2d defined by 

P = Y, P(vr)C T (vr) = E P (C T H) (8.29) 

and assume it irreducible and aperiodic. Then S holds and Theorem 8.6 applies with £ 
given by 

Xij = -^(rih) +^--^((^15(1)^) + ^-15(1)^)), (8.30) 

with S(l) the reduced resolvent of P at 1 and, for j 6 {1, . . . , d}, rj = ^2 reI± rj(r)|r) 6 C 2d . 
Proof: In this case, (8.15) reduces to 

(t' ® tt'\N r®i) = (t'ICV) t) P (tt'), (8.31) 
so that we can write with rji = I 71 ")) 

N T = ^p(vr)C T (7r)^| ??1 )(7r|. (8.32) 

7T 

Accordingly, for any £ G C 2d Cl n l, we have 

W T £ = C(0®lm>, with C(0 = Ef( 7r ) C ' T ( 7r )( 7r I^Cl«l' (vr|0cl«l e C 2d . (8.33) 

Hence, any eigenvector Vl' with eigenvalue e ld , 9 G R, needs to be of the form ^ = ^ <8> ??i 
with 

e*^ = p(^)C T (vr)^ = PV- (8-34) 

The matrix P being bi-stochastic, irreducible and aperiodic, there exists only one solution 
to (8.34), given by ip = ^2 re j ± \t) and e %e = 1, which shows that S holds. 

The expectation v and correlation matrix £ can be obtained from Theorem 6.6 in [21]. 
Indeed, under our assumptions, Lemma 8.3 shows that the process {Tj(uj)}j = \^.. tn is a 
Markov chain on I±, with kernel P = K p (C T (ui)) and initial distribution po( r o) = K 



2. 

TO I 



T(r„,. . . ,ti,t ) = |a ro | 2 E p(7Ti)p(7T2) ■ ■ ■ ,p(7Tn)('Tn|C(7r n )r„_i) ■ ■ ■ (ri|C(7ri)r ) 

7ri,...,7r n ef2 

= P(r n ,r n _i)---P(ri,r )po(ro). (8.35) 
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The aforementioned result provides the characteristics v and £ (8.30) of the functional 
central limit theorem for the Markov chain {Tj(UJ)}j = i r ^ jH corresponding to the random 
variable S n (uJ) = Y%=i r ( r j(^))- H 
Proof: With (3.57), (8.10) reads 

T(r n , . . . ,ti,t ) = |a T() | 2 ^2 P{^i)p(^2) • • • ,p(n n ) x (8.36) 

7ri,...,7r n eSl 

x(r n |C((7 E „-i r(T3) (7r n ))r„_i)-- - (ri|C((7o(7ri)))ro), 
where for all j > 1, thanks to the fact that a x is measure preserving, 

Setting P(t',t) = E p ({t' '\C r "(w)r)) and p (r) = |a r | 2 , we can write 

T(r n , . . . , n, r ) = P(r n , r n _i) • • • P(n, t )p (t ), (8.38) 

which proves the claim. M 

Remark 8.8 Actually, a strong law of large numbers holds in this case, i.e. liim^oo S n (U)/n 
r almost surely. 

9 Uncorrelated example 

In this last section, we briefly present two cases where the random coin matrices are chosen 
in an uncorrelated way, in order to complete the picture. In a sense, it can be viewed as 
the limiting case where the representation a of Z d is such that the periodicity lattice T is 
infinite. This is the complete opposite situation to the one considered in [21], where all coin 
matrices where identical in space, at all time step. Nevertheless, the methods developed in 
that paper apply here too. 

We recall some notations used in Section 2.1 in [21]: let x s = ^2j=\ r ( T j)i x 'j = 
Ylj=\ r ( r j)> then the generic term in Lemma 2.4 reads 



(r' s -i\C* s (x' s ) t')(t s \C s (x s ) r a _!) = (r>\C s (x> s ) t^X^C^) r s _x) (9.1) 

= (t s ®t' s \(C s {x s )®C s (x' s )) T^i^T^). (9.2) 

Let us introduce the unitary tensor product 

V s (x,y) = C s (x)®C s (y) in C 2d <g> C 2d . (9.3) 

Now consider the set of paths G n (K) in I? d from the origin to K = (J^j e ^ 2d v ^ a tne 
(extended) jump function defined by 



R:ll^, R['>\ = l*j>), (9.4) 
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that is paths of the form (Ti, ■ ■ ■ , T„_i, T n ), where T s = ^JJ G I±, s = 1, 2, . . . , n, and 

E"=i R ( T s) = K. For s>2 let X s = Y%Z\ R(Tj), while X x = 0. This last condition states 
that we start the walk at the origin. 

With these notations, we consider the complex weight of n-step paths in T? d from the 
origin to K, with last step T, defined by 

W£(n)= (T\V n (X n )T n ^)---(T 2 \V 2 (X 2 )T 1 )(T l V 1 (0) X o), (9.5) 

(T 1 ,...,T n _ 1 )ei% n ~ 1 B .t. 

with xo defined via the decomposition 

yo = «t|t) => Xo = ( A)®W = a ra T '|r r'). (9.6) 

re/± (r,T')e/| 

The expectation of this complex weight is the key quantity to analyze the averaged char- 
acteristic function (3.13), see [21]. Under certain assumptions on the distributions of the 
matrices Cj(x) G Q, fl finite, for simplicity, some cases can be readily studied using this 
method. 

Assumption A: 

a) The matrices V^(X) are distributed so that 

n 

nV^(X n ) = Z n ,C-i(*n-i) = Zn-ir-- ,V?(X ) = Z ) = P(yj"(^) = Zj). (9.7) 

b) The expectation E(V]?(X)) is independent of the position X: 

Qk= Y ZF{V?{X) = Z) = EW{X)). 

Assumption A is clearly satisfied in the following cases: 

Case 1: Assuming that the distributions of the matrices C s (x) are i.i.d in time and 
position, requirement a) is satisfied with P(V^(X) = Z) independent of j. Moreover, 
P(V w (x,y) = Z) = F (Z), for all x ^ y, and P(y w (x,x) = Z) = F D (Z), for all x. Further 
assuming 

Y ZF (Z)= Y ZF D (Z) = Q, (9.8) 

we meet requirement b) as well. 

Case 2: The following holds: 

i) For X € 7? d , V^(X) is a Markov chain in time on Q (g> Jl with initial distribution px 
and transition matrix Px- While for X ^ Y, the random variables V^(X), V UJ (Y) are 
independent. 

ii) The jump function R : l\ — > I? d is one to one and any X G I? d can only be reached at 
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most once on {£ re/ 2 ari? ( T )' "r G N} C l? d along any path X s = Y^ s j= \ R(Tj), s G N. 
hi) E(y/(A)) = Zzemn Zipx^Z) = Qj is independent of X, for any j G N. 



Under assumption A, we get the following expression for the expectation of W^(n) 



n-l 



S(w£(n))= ^ <T|Q n T n _i) JJ^ig^-iXTilQixo)- (9-9) 



(Ti,-,T n _ 1 )e/2™- 1 s.t. J=2 



'1.- .-<n-lJ*= J ± 
(T 1 ,..,T„_ 1 ,T)6G„(/f) 



Now we proceed as in [21]. Introduce the vectors in C 4 ^ 2 ~ C 2d <8> C 2d with Y G T 2d and 
n > 

*«on = E E e< ™ r w£wi r > and *o= E^it). (9.io) 

Using the notation 

D(y)= e iyR ( T ) |T)(T|, with F G T 2d and M k (Y) = D(Y)Q k , (9.11) 
re/2 

we obtain the following expression for the expectation 

E(* n (y)) = M n (y)M n _ 1 (y)---M 1 (y)* - (9.12) 

We get the following expression for the expectation of characteristic function (Proposition 
2.9 in [21]) 

Wi=/ (^^(y^M^^y^-.-M^y)^)^, (9.13) 

where 



*i = E i T > = E i r ^ r > and Y " = ( y v v ) € 

T&H ± rel± 



v>2d 



(9.14) 



At this stage the exact dependence of the matrix Mj on time j becomes crucial. In 
Case 1, Mj = M for all j, so that we are directly lead to the asymptotic study of 

f (*i|M"(y)* )^, (9.15) 

as in [21], which allows us to get diffusion properties and deviation estimates as in sections 
4, 5, 6, provided M(Y V ) satisfies the required spectral properties. 

In order to deal with Case 2 for non stationary initial distribution px, an analysis of 
the large j behavior of Qj based on the spectral properties of the transition matrix Px is 
in order. This should provide the necessary information to reach similar conclusions as in 
Case 1. 
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